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Recent cantilever magnetometry measurements of annular micron-size samples of S^RuCU have 
revealed evidence for the existence of half-quantum vortices (HQVs) in this material [l]. We pro- 
pose to look for HQVs in transport measurements and calculate magnetoresistance of a perforated 
superconducting film close to the transition temperature in the presence of HQVs. We analyze the 
dependence of magnetoresistance on the thermodynamic stability of HQVs which according to [T] 
can be varied by the application of an in-plane magnetic field and point out features which may 
help to identify them. 



Introduction. Half-quantum vortex (HQV) is a topo- 
logical defect in a two-component superfluid system 
which is defined as a state in which one of the components 
has an extra quantum of vorticity relative to the other. 
3 He-A, where the role of the two components is played by 
the spin projections, has been a long standing candidate 
in search for HQVs [3] however up to date the convinc- 
ing evidence for the existence of HQVs in that system 
is lacking [3j. More recently cantilever magnetometry 
measurements in micron-size Sr2Ru04 (SRO) samples, 
a candidate spin-triplet superconductor with the pairing 
symmetry which is believed to be analogous to that of 
3 He-A, have revealed evidence for the existence of HQVs 
in this material [T]. 

Recently revived interest in HQVs stems from their po- 
tential to host Majorana modes which in turn can provide 
a platform for topological quantum computing |9j . From 
a practical point of view the ability of SRO to support 
HQVs, if established, would provide a definitive proof 
for the unconventional nature of its pairing state. It is 
thus desirable, for both reasons, to have more experi- 
ments which could test the existence of HQVs in SRO. 
One possibility is to use scanning flux imaging employing 
either Hall [HIE] or SQUID [5J[TU] probes. However these 
techniques are technologically quite demanding and time 
consuming and not as readily available in experimental 
labs as e.g. transport probes. 

While most of the experimental work on SRO up to 




FIG. 1: (a) Proposed setup for detection of HQVs. Relevant 
dimensions are discussed in the text, (b) Elementary block 
model used for the calculation of magnetoresistance. See text 
for details. 



date has focused on bulk samples, recent technological 
advances made it possible to grow c-axis oriented epitax- 
ial films of SRO with thickness comparable to the zero- 
temperature coherence length and a reasonable transition 
temperature It is known that in superconducting 

films with nontrivial topology such as perforated films 
visualized as a collection of regularly spaced openings 
in an otherwise uniform medium ("swiss cheese" geom- 
etry, see Fig. [l^), variation of external magnetic field 
can produce resistance oscillations in the vicinity of the 
transition temperature [T21I15] . It is often a case that 
the periodicity of the most pronounced resistance oscil- 
lations corresponds to a flux quantum of the applied field 
through the smallest opening suggesting that such oscil- 
lations are caused by vortex entrances. Here we consider 
qualitative features of magnetoresistance (MR) of perfo- 
rated superconducting films in the presence of HQVs and 
point out to the signatures of HQVs in MR curves. 

Use of thin films in search for HQVs has a distinct ad- 
vantage that screening effects which in the bulk hinder 
the thermodynamic stability of HQV relative to that of 
a full vortex will be substantially reduced. The perfo- 
rated geometry which can be viewed as a collection of 
many small openings also offers advantage of relatively 
straightforward contact attachment which can be quite 
challenging for a small isolated SRO annulus. The down- 
side of a resistive measurement is that it should be done 
at temperatures close to T c where the phase stiffness of 
a spin current approaches that of a charge current [16] 
which makes HQVs less stable. However one of the key 
findings of ref. [T] is that the stability region of features 
attributed to HQVs can be expanded by applying in- 
plane magnetic field. Coupling to the in-plane field can 
provide a mechanism for stabilization of HQVs and hence 
can be used to both induce and identify their presence in 
the MR of thin films. 

Magnetoresistance of perforated superconducting films. 
Experimentally MR of thin perforated superconducting 
films shows an oscillating behavior with a monotonically 
decaying envelop caused by a global suppression of the 
superfluid density; the most pronounced oscillation pe- 
riod corresponds to a flux quantum of the applied field 
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through a smallest opening [51 [T^HHl IT7] . In high-quality 
arrays a fine structure such as less pronounced additional 
MR minima in each period has been observed and corre- 
spond to a flux quantum shared by several lattice cells. 
Our task is to show how to distinguish between these 
regular features and the ones induced by HQVs. 

A necessary condition for a film to produce observable 
MR oscillations is that it should be in a resistive state. 
One of the factors which controls the temperature range 
of the resistive state is the thickness of the film. While 
we do not see fundamental limitations on the film thick- 
ness apart from those set by obvious superconducting 
length scales, use of thinner films will broad the resistive 
transition thus facilitating the observation of MR oscil- 
lations. We should also point out that well defined MR 
features can be observed only in periodic arrays. After 
analyzing available data from different groups authors 
of ref. [14] concluded that one of the conditions for the 
observation of MR oscillations is that the width of the 
constriction separating neighboring openings should be 
less than 4£(T). However from the point of view of ther- 
modynamic stability of HQVs for films with thickness 
w < A it is beneficial to reduce the constriction width 
below the Pearl penetration depth A 2 /u> which for SRO 
where A/£«2.5 can be larger or smaller than 4£(T) de- 
pending on the film thickness. Thus a priori to improve 
the possibility of seeing HQVs in perforated SRO films 
for films with w < A(T) openings can be spaced further 
apart (but not more than 4£(T) for MR) while for thicker 
films openings should be spaced as close as possible. 

Theoretical model. We now proceed to the calculation 
of magnetoresistance (MR) of a perforated superconduct- 
ing film close to T c . Existing approaches to this problem 
fall into two groups. In the first group one considers a 
periodic array of openings and calculates variation of T c 
induced by the applied field H; it is then argued that the 
MR will trace the dependence T C (H). Such approach has 
been adopted e.g. for the proximity induced Josephson 
link arrays 12J, for a regular square network of supercon- 
ducting aluminum [13] and for thin amorphous TiN films 
which show superconductor- insulator transition [14] , 

Approaches in the second group are based on the ob- 
servation that the most pronounced features in the MR 
are correlated with a quantized value of the flux though 
a single opening so that instead of considering the resis- 
tance of the whole lattice one can focus on determining 
the resistance of a much smaller ("elementary") block 
which usually consists of just one [17] or few [12] lattice 
cells. Given a particular block structure it is then neces- 
sary to establish a dissipative mechanism responsible for 
a finite resistance. Close to T c where £ becomes compara- 
ble to the width of constriction between the openings the 
mechanism responsible for the oscillatory part of MR is 
provided by phase slips generated by thermally activated 
vortex jumps in and out of openings. It is then argued 
that such phase slips are described by the Ambegaokar- 



Halperin (AH) model of a thermally activated resistance 
of a heavily damped Josephson junction [18J. 

AH considered RCSJ model of a Josephson junction 
supplemented by a thermal voltage noise. They showed 
that this problem is equivalent to the problem of Brow- 
nian motion of a particle in the "washboard" potential 
and, in the case of large damping, obtained I-V curves 
of the junction. In the limit of vanishingly small currents 
the ratio of the junction's resistivity to its normal state 
value R n is given by 

R/R n = Io 2 (AE/2k B T) (1) 

where Iq is a modified Bessel function and AE is the 
activation barrier. 

The AH model has been successfully applied to a va- 
riety of physical situations seemingly beyond its original 
scope. In particular, assuming that vortex passages can 
be treated as phase slips and choosing appropriate acti- 
vation barrier Tinkham has applied eqn. ([I]) to explain 
the width and the shape of the irreversibility line in some 
cuprate superconductors [T5]. More recently Sochnikov 
et al. 17 have successfully modeled magnetoresistance 
of nanopatterned superconducting LSCO arrays using 
eqn. ([!]) with the activation barrier determined by the en- 
ergy level spacings between vortex states of array's open- 
ings. We will adopt a similar approach and use the AH 
model to point out features in the oscillatory part of MR 
which are brought about by the presence of HQVs. 

One can argue that application of the AH model to 
multiply connected geometries in the magnetic field re- 
quires further justification since eqn. |l]) was derived un- 
der the assumption of vanishingly smalls currents while in 
the magnetic field, even in the absence of transport cur- 
rent, there generally exists non-zero screening current. 
To justify the use of eqn. ([I]) for that case we consider a 
pair of Josephson junctions threaded by flux $ (dc-squid) 
with each junction described by the RCSJ model with 
thermal voltage noise. By writing down a set of coupled 
equations for phase drops 6\ and 9 2 of the junctions it 
can be shown that under the assumption of negligible flux 
noise this model describes a single junction with phase 
drop (02 — 0i)/2 and a flux-dependent activation barrier 
and is thus equivalent to the AH model. 

Thus we have reduced the problem of calculating of 
MR of a perforated superconducting film to a calcula- 
tion of a flux-dependent activation barrier of an elemen- 
tary block. Further simplification is reached by noticing 
that provided the ratio of the Josephson energy to the 
superfluid density is not very small the shape of the ac- 
tivation barrier of a dc-squid is well represented by the 
spacing between energy states with different vorticity of 
a junction-free loop which can be checked through their 
graphical comparison. This observation allows us to re- 
place a dc-squid in the elementary block with a junction- 
free loop with uniform phase winding. 

For the reasons which will become clear later we con- 
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FIG. 2: Magnetoresistance of a perforated film around 
^/$o = 1/2 in the presence of HQ Vs for different values of the 
in-plane field. The curves are offset for clarity and correspond 
to the in-plane field H x determined by, from the bottom to 
the top, h(H m )/eo = 0.04, 0.05, 0.07, 0.09 and 0.10. All curves 
are evaluated at 7 = 4/5 and eid/^o = 0.05 for which H x q is 
given by h(H x0 )/e = 0.05. 

sider an elementary block which is composed of two lat- 
tice cells and is represented by two overlapping loops with 
circumference L each, the overlap length £ and cross- 
section A as on Fig. [T|d. In the limit T — » T c the ra- 
dial variations of the order parameter can be neglected 
and each loop which can support HQVs or, more pre- 
cisely, half-integer fluxoids should be characterized by 
two winding numbers which we denote as n s j and n sp i, 
i = 1,2 (20] [21]. In this notation full vortices are de- 
scribed by n si G integer and n spi = while for HQVs 
n si G half-integer and n sp i — ±1/2. 

We are now in the position to consider energy levels of 
an elementary block which, according to the discussion 
above, can be used to compute activation energy barriers 
needed for the calculation of the oscillating part of MR. 
We notice that close to T c the screening is weak so that 
instead of Gibbs potential one can use free energy; the 
reduced screening also allows us to ignore the inductive 
coupling between the loops. Denoting winding numbers 
of each loop n s ,, n sp i and introducing n s ± = (n s i±n S 2)/2 
and, analogously, n sp ± = (n sp i ± n sp 2)/2 the free energy 
of the two-cell elementary block with each loop pierced 
by applied flux <£> is given by 
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where e = A^/(4tt\ 2 (L - t)) and 7 = (L-l)/(L+t) < 
1 is a numerical factor determined by the geometry of the 
elementary block. 

The free energy specified by eqn. ^ contains two con- 
tributions. The first contribution proportional to the 
term in the brackets - the kinematic term - is the ki- 
netic energy of charge and spin currents which depends 
on p s and p sp - superfluid and spin superfluid densities, 
which characterize stiffness of charge and spin currents 
and satisfy the inequality p sp / p s < 1 |16| . If the screening 



is ignored the latter condition implies that the kinematic 
term alone can provide the stability of HQV [20] [21]. 
However upon approaching the transition temperature 
Psp I ' Ps ~ ► 1 so that the kinematic contribution by itself 
cannot guarantee the thermodynamic stability of HQV. 

The second contribution in eqn. ^ denoted as Ae HQV 
contains terms which are explicitly flux-independent 
e.g. those coming from spin-orbit interactions. Recent 
cantilever magnetometry measurements of micron-size 
SRO annuli [1 , if interpreted in terms of HQVs, imply 
that the magnitude and the sign of this term can be con- 
trolled by the in-plane field - the magnetic field directed 
along the ab-plane which in our geometry coincides with 
the plane of the film and hence does not contribute to 
the flux $. While the origin of such in-plane coupling is 
still unclear, it may be attributed to the effect of kine- 
matic spin polarization according to which the HQV state 
posses additional spin polarization not present in the full 
vortex state [21 j - Coupling of the kinematic spin po- 
larization to the in-plane field requires that the pairing 
symmetry of SRO is such that the axis of "easy" spin 
polarization - so-called equal-spin-pairing axis - has a 
non-zero in-plane component. 

Generalizing the in-plane field coupling model pro- 
posed in [1] for a two-loop elementary block Ae HQV can 
be written as 



Ae H 



eid(\n sp i\ + \n sp i\) - h(H x )(n spl + n sp2 ) (3) 



where €m > is a term which originates from the spin- 
orbit energy difference between the half-quantum and full 
vortex states and h(H x ) is a monotonic function of the 
in-plane field H x ; notice that the coefficient in front of 
h(H x ) is proportional to the sum of spin currents in each 
loop in accordance with the kinematic spin polarization 
model [21 . As observed in pQ for micron-size SRO annuli, 
contribution Ae H Qv which is non-zero only in the HQV 
states, becomes negative as the in-plane field reaches 
a sample-dependent value H xQ which lies in the range 
10G-200G thus inducing the thermodynamic stability 
of HQV. In our model H x q is defined by the condition 
h{H x0 ) — eid- 

The full set of energy levels described by eqn. ^ is 
periodic in flux and is degenerate with respect to the in- 
terchange of the two cells in the elementary block; we 
will ignore this degeneracy since in real-life situations it 
will be removed by the presence of the sample's bound- 
ary. Complexity of the energy levels grows substantially 
with their energy and depends on the value of param- 
eters entering eqn. |2]). In the absence of the in-plane 
field the lowest energy level consists of three full vortex 
states. Denoting each state by a set of four winding num- 
bers {n s i,n sp i,n s i,n sp2 } and setting p sp /p s = 1 as ap- 
propriate for temperatures close to T c , the lowest energy 
level is realized by {0,0,0,0} for $/$ G [0, (1 + 7)/4], 
{1, 0, 0, 0} for $/$ G [(1 + 7 )/4, (3 - 7 )/4] and {1, 0, 1, 0} 
for $/$q G [(3 — 7 )/4, 1] which corresponds to no vortices, 
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one vortex and two full vortices in the elementary block. 

As the in-plane field reaches value H x0 defined above 
the lowest energy level acquires contributions from states 
{1/2,1/2,0,0} and {1/2,1/2,1,0} which correspond to 
an HQV in one of the cells. Notice that states with 
HQVs in each cell will either not couple to the in-plane 
field (state {1/2, 1/2, 1/2, -1/2}) or will contribute to the 
ground state level only starting from a higher value H x \ 
of the in-plane field (state {1/2,1/2,1/2,1/2}) defined 
by h(H xl ) = e ld + e„(l - 7 )/4. 

Results. The oscillating part of magnetoresistance is 
calculated using eqn. ([lj choosing the activation barrier 
as a spacing between two lowest energy levels of model 
([2]); the result is shown on Fig. [2] for different values 
of the in-plane field. Zooming in the neighborhood of 
$/$o = 1/2 we see that for in-plane fields H x < H x o 
magnetoresistance has a dip at <&/$o = 1/2; this dip 
is produced by a configuration with one full vortex per 
two-cell elementary block and can be observed in con- 
ventional superconducting films, see e.g. [12) . Upon in- 
creasing the in-plane field above H x q states with one 
HQV per block produce satellite dips located symmet- 
rically around $/$o = 1/2 while the central minimum 
becomes shallow and eventually disappear. Unlike the 
minimum at $/$o = 1/2 such behavior is characteris- 
tic only to the films which are able to support HQVs. 
For in-plane fields just above H x o location of the satel- 
lite HQV minima $ HQV is field-independent and is given 
by |$ HQ v/$o - 1/2| = (1 - 7)/4- 

The presence of higher energy levels can be accounted 
for by using a procedure suggested by Sochnikov et 
al. [T7]- This procedure consists of replacing AE cal- 
culated on the two lowest energy levels with the Gibbs 
average (AE) calculated using full spectrum described by 
([2]). We have checked that such procedure while able to 
smooth out magnetoresistance features does not change 
their evolution with the in-plane field as described above. 

It should be noted that conventional vortex configura- 
tions can also produce magnetoresistance minima which 
are located at the rational values of the applied flux. 
These features, which are usually less pronounced, cor- 
respond to full vortex states of larger elementary blocks 
and can only be observed in high-quality films [131 HI • 
However, unlike minima produced by HQVs, minima due 
to conventional vortices will not evolve with the in-plane 
field in the way described above. Hence the in-plane field 
evolution of HQVs as observed in [T] can be used to iden- 
tify their magnetoresistance signatures. 

Conclusion. We have proposed to search for half- 
quantum vortices in transport properties of perforated 
superconducting films in the presence of magnetic field. 
In particular, we have calculated magnetoresistance of a 
perforated superconducting film assuming that the film's 
material is able to support half-quantum vortices as has 
been recently experimentally suggested for Sr2Ru04 in 
cantilever magnetometry measurements pQ. We have 



shown that the oscillating part of the magnetoresistance 
should have well defined signatures of HQVs which, given 
the findings of ref. pQ can be controlled by applying the 
in-plane magnetic field. 
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